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Design  Factors  in  the  Tuned 
Synthetic  Aperture  Radar 

I.  INTRODUCTION 

Since  the  late  1930s,  one  of  the  most  impressive  areas  of  technological 
development  has  been  that  of  radar.  However,  in  spite  of  our  present  expertise, 
the  operational  capability  to  detect  and  track  slowly  moving  targets  accurately, 
on  or  near  the  ground,  from  a swiftly  moving  airborne  platform  does  not  exist. 

The  Tuned  Synthetic  Aperture  Radar  (TUSAR)  technique  aims  at  achieving  that 
capability.  The  development  of  TUSAR  would  have  its  application  in  long-range 
battlefield  surveillance.  If  mounted  in  a C-130,  for  instance,  such  a radar  could 
inform  a ground  commander  of  the  entire  situation  of  moving  men,  armaments, 
and  supplies,  and  never  get  closer  than  40  miles  from  the  action. 

The  mathematical  and  theoretical  foundations  for  the  Tuned  Synthetic  Aperture 
Radar  have  been  laid. 1 One  realization  of  the  Tuned  Synthetic  Aperture  Radar 
processing  is  evaluated  in  this  paper.  * The  evaluation  requires  a computer 
program  of  sizeable  magnitude.  The  program  can  be  broken  down  readily  into 
several  sub-programs.  The  sub-program  described  herein  determines  the 


(Received  for  publication  24  August  1977) 

1.  Schindler,  J.K.  and  Coggins,  W.  B.  An  Airborne  Radar  Technique  for  Moving 
Target  Detection,  Location  and  Tracking.  AFCRU-TR-73-0719. 

* This  work  was  conducted  while  one  of  the  authors  (RAL)  was  participating  in 
a cooperative  research  program  from  the  US  Air  Force  Academy. 
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number  of  radar  outputs,  particular  outputs,  and  number  of  data  sets  required  for 
processing.  The  development  of  this  sub-program  came  about  in  three  stages, 
the  first  of  which  was  determining  the  exact  function  to  be  performed  and  the  use 
to  which  the  results  will  be  applied.  The  second  stage  involved  developing  pro- 
gram logic,  and  the  third  stage  involved  using  actual  data  in  the  program  to  deter- 
mine trends  which  might  be  useful  in  later  design  stages. 


1.  I’ltCM.K  \M 

The  digital  AMTI  processing  employed  by  the  Tl’SAR  technique  is  to  be  evalu- 
ated with  a truck-borne  test-bed  radar  operating  at  X band.  The  inphase  and 
quadrature  components  of  the  received  signal  are  digitized  and  processed  in  a sig- 
rial  processing  oriented  MAP300  mini-computer  for  doppler  frequency  selection. 
The  signal  is  passed  through  a Fast  Fourier  Transform  (FFT)  of  as  yet  undeter- 
mined size,  and  successive,  contiguous  doppler  frequency  outputs  are  linear!-,  and 
coherently  combined  and  detected  to  yield  target  detection  data.  Specifically,  as 
tiie  radar  test  vehicle  moves,  the  complex  amplitude  of  the  received  signal  from 
each  transmitted  pulse  is  observed.  When  the  returns  from  N pulses  are  received, 
they  are  batch  processed  by  an  FFT,  yielding  N doppler  frequencies  at  the  output. 
This  process  is  repeated  over  M successive  but  not  necessarily  contiguous  time 
periods  of  N pulses,  allowing  one  to  track  the  target  by  following  its  linear  doppler 
frequency  modulation  over  the  N x M pulse  period  of  time. 

In  general,  the  processor  has  NxM  doppler  frequency  data  points  available. 

If  only  slowly  moving  targets  are  considered,  a small  portion  of  these  N x M data 
points  are  needed.  For  instance,  selecting  the  initial  doppler  frequency  of  the 
target  in  the  range  - . l i'cm  to  4 . lu>cm  (w  m = maximum  doppler  radian  frequency 
of  the  ground  clutter)  and  a maximum  frequency  change  of  . 2tu  over  the  total 
radar  integration  time  means  that  moving  targets  must  appear  in  a doppler  fre- 
quency band  of  width  . 4t  . The  Nyquist  sampling  theorem  requires  a pulse 
repetition  frequency  of  at  least  twice  the  highest  frequency.  That  means  the 

available  FFT  outputs  cover  a bandwidth  of  2iu  or  more.  Needing  only  . 4di 

cm  ° J cm 

bandwidth  to  observe  slowly  moving  targets  allows  us  to  ignore  at  least  80%  of  the 
possible  FFT  outputs.  The  resulting  simplification  in  the  processing,  through 
pruning  of  the  FFT  algorithm,  will  be  utilized  here. 


*The  MAP300  is  a high-speed  digital  signal  processor  manufactured  by  CSPI, 
Burlington,  MA. 


i 


i 

I 


6 


1.  I II*  I INK  \K  lM>m.KR  KKK.M  KM  H VITKOMMVI'ION 


I 


i 


■« » . ..  ~ 


fho  tuned  Synthetic  Aperture  Kadar  technique  is  based  on  the  idea  that  within 
a limited  segment  of  the  along  track  path  of  the  radar,  the  doppler  frequency  shift 
of  both  constant  velocity  moving  targets  and  stationary  clutter  is  linear  (figure  1) 


figure  1.  A Typical  Doppler  History,  Showing 
Region  of  l inearity 


This  approximation  improves  as  the  doppler  frequency  bandwidth  of  the  moving 
target  and  the  limited  segment  of  along  track  path  decreases.  The  linear  fre- 
quency approximation  simplifies  the  process  of  detecting  a moving  target  to  pos- 
tulating the  slope  and  initial  frequency  of  the  target  and  then  incrementing  anil 
coherently  combining  the  doppler  frequency  output  of  the  successive  N point  1 1 1'.-, 
to  coherently  track  the  doppler  frequency  of  the  target.  This  is  the  baste  process 
assumed  in  the  development  of  a program  which,  given  the  frequency  of  entry  and 
rate  of  frequency  decreases  of  the  moving  target,  provides  a sequence  of  ffT 
outputs  that  coherently  tracks  the  moving  target  over  the  integration  time  of  the 
radar. 

t \n;  vm  k KKP.yi  km.ik.s 

As  indicated  in  figure  1,  the  instantaneous  doppler  frequency  of  a moving 
target  decreases  with  time.  Targets  will  then  exhibit  negative  frequencies. 

The  program  under  development  must  allow  for  this  problem  in  a system  where 
nominally  only  positive  frequencies  occur  at  the  fl'T  outputs.  The  sampling 
theorem  is  once  again  the  key.  figure  2 shows  the  sample  return  frequency 
spectrum  with  its  recurring  patterns  at  harmonics  of  the  radar  pulse  repetition 
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frequency,  i ■ Negative  frequencies  are  shifted  to  the  next  harmonic  (see 
Figure  3),  ami  read  as  positive  frequencies  just  below  r As  the  target's 
doppler  frequency  decreases,  it  moves  from  right  to  left  in  the  shaded  area 
commencing  at  A.  When  the  target  reaches  zero  doppler  frequency  it  jumps  to 
tprf,  then  continues  to  shift  from  right  to  left  as  though  passing  through  negative 
frequencies.  The  output  of  each  N point  FFT  is  N different  frequencies  ranging 
from  zero  to  i j.  Numbering  each  of  these  with  an  integer  value  from  1 to  N 
provides  a convenient  way  of  referencing  the  particular  outputs  to  bo  used,  as  well 
as  telling  us  how  many  of  the  outputs  are  needed  for  slow-moving  target  detection. 


Figure  2.  Received  Signal  Spectrum  With  Range 
of  Target  Frequencies 


Figure  3.  Received  Signal  Spectrum  With  Shifted 
Range  of  Target  Frequencies 
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Assuming  that  the  moving  targets  occupy  40'ru  of  the  maximum  elutter  frequency, 

i , and  that  clutter  is  barely  unambiguous,  (j  . ),  then  the  number  of 

cm  • B prf  cm 

KKT  outputs  required  are 


• 4j‘ 

cm  N 
“1prf 


. 4f  V 

— 1 L v 

c 1*111 


• ZN. 


Here  f is  the  radar  carrier  frequency,  V is  platform  velocity,  PRF  is  pulse 
repetition  frequency,  and  c is  speed  of  light. 


r*.  KRH.M  KM  A IMK  K\  \1.S 

The  frequency  and  time  scales  of  the  observations  are  conveniently  divided 
into  blocks  of  N/l’KF  time  units,  and  PRF/N  frequency  units.  These  units  are 
derived  from  the  fact  that  each  FFT  requires  N pulses,  and  each  interpulse  time 
is  l/l’HF,  Hence,  there  are  N'/l’HF  time  units  per  block.  Since  the  frequencies 
at  the  outputs  range  from  0 to  r and  there  are  N of  these  outputs,  there  must 
be  PRF/N  radians  per  second  per  block.  Figure  4 depicts  these  "L'ocks''  graphi- 
cally. They  are  intervals  of  interest  with  integer  assignments  and  co.i-'ant  width 
or  height. 


Figure  4.  Division  of  Time  and 
Frequency  Intervals  With  Postulated 
Target  Type  Superimposed 


These  blocks,  or  "bins,  " represent  the  possible  time-frequency  combinations 
that  a target  can  attain.  In  Figure  4,  a line  has  been  drawn  to  represent  a pos- 
tulated moving  target.  Kach  box  through  which  this  line  passes  indicates  an  FFT 
output  at  a certain  time  where  the  target  signal  will  appear.  Notice,  however. 


that  there  are  bins  which  are  barely  touched  by  the  target  and  which  should  per- 
haps be  ignored  since  the  contribution  to  improving  the  target  signal  to  clutter  of 
these  bins  is  minimal.  This  condition,  indicated  by  shaded  boxes,  means  that 
one  should  selectively  omit  certain  bins  from  the  output  processing.  This  pro- 
cess may  be  accomplished  by  a percentage  lop-off  procedure,  which  observes 
the  fraction  of  time  a target  spends  passing  through  a bin  and  ignores  those  bins 
that  contribute  little. 


f>.  LIMITS  ON  TUT.  I K \Cf 

The  basic  problem  at  hand  is  to  filter  clutter  and  improve  the  target  signal- 
to-clutter  ratio.  A decision  must  be  made  at  this  point.  Should  the  search  for 
a postulated  target  pass  through  a fixed  number  of  frequency  bins,  allowing  the 
number  of  time  intervals  to  change,  or  would  the  SCR  be  greater  by  tracking  each 
target  type  for  a given  number  of  time  periods,  allowing  the  number  of  frequency- 
bins  occupied  by  a target  to  vary?  Although  no  quantitative  analysis  was  per- 
formed, intuition  coupled  with  qualitative  arguments  suggested  using  a fixed 
number  of  frequency  bins  for  each  target,  regardless  of  the  time  frame  needed. 
Thus,  the  algorithm  was  designed  to  select  a star*  frequency  and  terminate  track 
after  a frequency  change  of  0.  &lcm  was  achieved. 


7 1*1  RPOSK.  OK  1'IIK  I’KOCR  \M 

The  purpose  of  this  program  is  twofold.  The  first  is  the  most  obvious;  it 
shows  the  exact  succession  of  FFT  outputs  needed  for  target  detection,  how  many 
are  required  for  each  time  period,  and  how  much  memory  is  required.  The 
other  purpose  is  more  subtle.  The  FFT  output  in  each  time-frequency  bin  is  to 
be  weighted  and  coherently  added  to  others  to  form  a single  detection  statistic. 

The  relative  weighting  factor  can  be  found  by  applying  signal -to-clutter  maximi- 

1 0 

zation  techniques  described  by  Coggins  and  Schindler.  ’ “ The  matrices  involved 
in  determining  the  maximum  signal -to-clutter  weighting  factors  are  dependent 
upon  the  sequence  of  FFT  outputs  which  are  assumed  to  be  combined.  The  pro- 
gram described  here  was  designed  to  evaluate  the  successive  FFT  outputs  re- 
quired for  processing.  Future  work  will  address  the  design  of  the  relative  weights 
required  for  coherent  processing  of  these  FFT  outputs. 


2.  Coggins,  W.  B.  and  Schindler,  J.K.  Processing  for  Maximum  Signal-to- 
Clutter  AMT1  Radars.  AFCRL-TR-74-0577. 
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How  are  specific  target  parameters  chosen  for  design  purposes?  Obviously, 
there  is  an  infinitely  large  set  of  unique  velocities  and  ranges  to  be  considered. 
Initially,  this  set  of  target  parameters  must  be  constrained  to  lie  within  certain 
bounds  determined  by  the  maximum  ground  speed  of  the  target  and  its  maximum 
and  minimum  range.  Then,  a finite,  discreet  set  of  target  speed  and  direction 
pairs,  as  well  as  a set  of  target  ranges,  must  be  selected  to  serve  as  design 
parameters.  Targets  whose  parameters  fall  in  the  gaps  between  the  selected 
design  parameters  must  be  detected  by  virtue  of  the  finite  resolution  of  the  pro- 
cessing. VVitli  the  speed  available  in  modern  digital  processing  equipments,  it 
appears  quite  possible  to  think  of  processing  200-300  general -target  types  in  the 
time  required  to  gather  312  pulses  of  raw  data  at  UHF  radar  frequencies. 


«.  PROGRAM  l.()C;i('. 

Having  provided  the  basis  for  the  program's  development,  it  is  time  to 
discuss  the  logic  involved.  The  program,  written  in  Fortran  (extended)  for  the 
CDC6600,  consists  of  five  basic  stages.  The  first  stage  is  a simple  interactive 
scheme  for  inputting  data  through  a remote  terminal  and  spotchecking  input  data. 
This  stage  includes  program  lines  100  through  330  (a  listing  of  the  computer  pro- 
gram is  given  in  Appendix  1.)  In  the  next  stage  (lines  3f!0  to  500),  the  radar  and 
target  parameters  discussed  by  Schindler  and  Goggins*  are  evaluated. 

l.ines  590  through  1280  contain  the  main  iteration  loop.  After  calculating  the 
doppler  frequencies  where  the  target  enters  and  leaves  a specific  time  interval, 
the  first  step  is  to  determine  whether  the  target  contributes  significantly  to  the 
output  power  from  an  FFT  filter  during  that  interval.  The  procedure  is  to  com- 
pare the  portion  of  the  frequency  axis  intercepted  by  the  target  with  the  total 
length  of  the  frequency  interval,  that  is,  the  bandwidth  of  the  FFT  filter.  If  a 
specified  percentage  (LOP)  is  not  met  or  exceeded,  that  frequency-time  interval 
is  eliminated.  Figure  5 illustrates  this  process  at  entering  and  exiting  sides  for 
slopes  for  less  than  -1.  This  case  applies  when  the  target  doppler  bandwidth 
during  the  time  observation  interval  N/P11F  exceeds  the  FFT  filter  cycle  band- 
width, PRF/N.  When  the  slope  is  greater  than  -1  (Figure  G)  the  time  excursion 
of  a target  is  compared  to  the  whole  time  interval,  N/PRF,  to  determine  accept- 
able targets.  In  each  figure,  line  segment  1 as  a percentage  of  line  segment  2 is 
compared  to  the  parameter  LOP  to  determine  acceptable  targets. 

Having  decided  what  bins  to  ignore,  the  program  goes  on  to  store  acceptable 
bin  numbers  for  later  output.  There  are  three  possible  cases  here.  Either  the 
time  interval  contains  all  positive  frequencies,  some  positive  and  some  negative 
frequencies  or  all  negative  frequencies.  After  recording  all  the  bin  numbers  for 
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Figure  5.  I .OP  Hate;  Slope<-l  Figure  6.  LOP  Hate;  -l<Slope<0 


the  interval  in  one  of  the  three  categories,  the  program  steps  to  the  next  time 
interval  and  repeats.  Lines  after  1 25)0  then  output  the  computed  information  in 
three  forms:  input  parameters,  specific  outputs  numbered  sequentially  and  within 
each  FKT,  and  number  of  outputs,  chops,  and  time  intervals  involved. 

<>.  INMCVmK  TKKNDS 


Once  completely  debugged,  the  program  was  ready  for  use  in  theoretical  case 
studies.  One  of  the  first  studies  performed  was  to  change  the  FFT  size  and  LOP 
rate,  holding  all  radar  and  target  parameters  constant,  and  observe  the  number 
of  FFT  outputs  required.  The  results  of  this  study  are  plotted  on  two  graphs. 
Figures  7 and  8.  Figure  7 shows  a linear  decrease  in  the  number  of  FFT  outputs 
required  as  a function  of  the  LOP  rate.  This  fact  indicates  that  we  can  freely 
change  the  LOP  rate  to  achieve  the  desired  number  of  outputs.  Figure  8 points 
out  two  interesting  facts.  First,  as  a function  of  FFT  size  and  hence  FFT  band- 
width, there  is  a point  where  the  number  of  bins  accepted  is  at  a minimum.  The 
minimum  occurs  when  the  average  frequency  change  during  an  FFT  observation 
time,  N/PRF  equals  the  FFT  bandwidth  PKF/N.  The  curve  is  parabolic  when  the 
abscissa  is  defined  in  log^  units.  The  second  interesting  fact  is  that,  at  the  min- 
imum, the  number  accepted  is  a sensitive  function  of  the  LOP  rate.  This  is 
manifested  in  the  spread  between  the  curves  at  the  minima,  as  compared  to  else- 
where on  the  curve. 

Another  interesting  effect  was  noticed.  As  the  FFT  size  increases,  the  num- 
ber of  time  periods  required  decreases.  The  overall  effect  is  to  produce  a 
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reasonably  constant  ti>tal  number  of  outputs  to  Lm*  stored  for  processing.  Should 
a designer,  for  reasons  of  resolution  or  computation  speed,  want  to  change  the 
!•' l-’T  size,  he  can  he  assured  such  a change  will  not  drastically  affect  the  mem- 
ory allocations. 

On  the  other  hand,  additional  studies  indicated  that  there  can  be  as  much  as  a 
70":  change  in  the  memory  requirements  when  target  and  platform  velocities  are 
adjusted.  This  means  that  the  ranges  of  all  parameters  will  necessarily  have  to 
be  decided  upon  just  before  the  memory  structure  is  designed. 

Also  of  interest  was  the  fact  that,  in  the  0"',.  1.01*  case,  changing  the  FFT  size 
produces  "glitches"  in  the  parabolic  curves  of  Figure  8— areas  of  spikes  where  the 
number  of  outputs  accepted  suddenly  increases.  This  is  due  to  the  digitization 
process  which  tries  to  fit  the  straight-line  doppler  history  into  a series  of  inter- 
vals. These  spikes,  however,  tend  to  be  smoothed  at  higher  l.Ol’  rates. 
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10.  SI.OPK.  m'  Sl/K  CAMT  RATIONS 

As  shown  in  Figure  8,  there  is  an  FFT  size  where  the  fewest  number  of  out- 
puts ure  accepted.  This  occurs  when  the  average  slope  of  the  target  doppler 
history  is  negative  unity.  The  doppler  rate  or  slope  is  proportional  to  the  change 
in  target  doppler  frequency,  Af.  observed  in  a time  period  compared  to  the  size  of 
the  FFT  frequency  interval,  FI.  That  is,  the  slope,  m,  is  given  by 

m = Af/ FI  . 

Now  Af  can  lie  expressed  as  the  target  doppler  rate  times  length  of  time  considered. 
The  frequency  interval,  FI,  is  merely  PRF/N.  Thus 


m = — — — / PR  F/  N = — -N“^ 
2rr  PRF  2tt(PRF) 


A reasonable  set  of  target  doppler  rate  parameters  lies  between  1.21ti>  and 

cm 

. 81;i>  , where  i>  is  tlie  maximum  doppler  rate  of  stationary  clutter  elements. 

This  range  corresponds  to  ground  targets  traveling  at  a speed  which  is,  at  most, 

10"', i of  the  radar  platform  speed.  It  is  thus  reasonable  to  use  ib  for  the  target 

cm  b 

doppler  rate  in  the  above  formula  form,  since  i is  at  the  approximate  mean  of 
the  target  doppler  rates. 

From  (1) 


where  R{  = range  to  the  target. 
Then 


4nf  V2N2  .. 

m — 2tt(PRF)“ 


o ■> 

2f  v N‘ 


RtC  (PRF)* 


- 


wrr,,t  1 


To  find  the  relative  minimum  in  the  curve  of  accepted  outputs  versus  I’ FT 
size  (Figure  8),  allow  m = -l  and  solve  for  N: 


The  value  of  N specified  here  gives  the  optimal  match  of  FFT  bandwidth  of  ground 
clutter  doppler  excursion  during  the  N/PRF  integration  time.  Individualities  in 
the  moving  targets  doppler  properties  may  change  slightly  where  this  minimum 
occurs. 

II.  PROCESSOR  KKyiURKMKNTS 

Numerous  calculations  and  estimates  were  made  to  insure  that  the  entire 
processor  function  could  be  performed  within  the  time  and  memory  limits  of  the 
truckhorne  CSP-30/  MAP300  processor,  and  would  later  be  practical  in  an  air- 
craft. To  start  with,  the  processing  time  requirements  in  the  truckhorne  case 
were  examined.  There  are  basically  four  functions  for  the  central  processor  to 
perform:  multiplication,  addition,  memory  retrieval  and  counter  update.  The 
maximum  time  allowable  for  any  processing  of  a given  data  set  is  going  to  be  the 
time  it  takes  to  collect  N data  points.  After  calculating  the  computational  times, 
a good  estimate  would  be  to  double  that  time  to  allow  for  program  execution.  The 
chart  given  in  Figure  9 shows  these  times.  The  times  are  based  upon  the  truck 
case  requiring  summation  of  50  complex  numbers,  each  multiplied  by  its  own 
weighting  factors.  It  is  assumed  that  the  radar  PRF  is  1831  Hz,  and  512  point 
FFTs  are  required  in  a test  for  200  possible  target  doppler/ doppler  rate  com- 
binations. The  reader  should  note  that  these  are  not  intended  as  exact  numbers, 
but  rather  are  approximations  to  prove  that  processing  time  will  not  be  a limiting 
factor.  The  FFT  process  will  be  run  simultaneously  on  a different  part  of  the 
processor,  requiring  only  4.  5 ms  for  a 1024  complex  point  FFT,  and  will  be  idle 
98.4%  of  the  time.  Overall,  the  speed  with  which  the  computations  take  place 
allows  a vast  number  of  target  possibilities  to  be  examined  without  stressing  the 
real-time  processing  capabilities  of  the  CSP-30. 

The  other  major  constraint  on  computational  abilities  is  the  size  of  the  mem- 
ory. The  CSP-30  has  an  existing  32K  word  memory  available,  each  word  having 
1G  bits.  Figure  10  shows  a breakdown  of  the  various  memory  requirements. 
These  numbers  assume  the  truckhorne  case  of  40  time  periods,  200  tasks,  50 
summations  per  target,  and  only  l memory  access  value  needed  for  each  data 
point.  They  also  assume  a rectangular  matrix  structure  for  memory,  for  which 
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FUNCTION 

CSP  30 
UNIT  TIME 

TOTAL  TIME 

- zoo  multiplies 

1 100  nice 

220  |.uc 

i 100  AOOS 

400  ni«c 

40  ,.»«c 

100  ME  MOST  RETRIEVALS 

300  n*«c 

30  fittc 

50  COUNTER  UPOATES 

200  nuc 

10  F»«c 

Figure  9.  Computation  Times 

300  it** c 

200  REPETITIONS 
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DATA  COLLECTION  TIME  (312  pt  FFT)  280  niMC 
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OAT  A POINTS 

40 

TIME  PERIODS 

103 

POINTS  SAVEO  EACH  FFT 

* 2 

WORDS  PER  COMPLEX  WORD 

6240 

TOTAL  WOROS 

MEMORY  REFERENCE  SEQUENCES 

50 

POINTS/SUMMATION 

« 200 

SUMMATIONS 

10,000 

TOTAL  WOROS 

woros 


OATA  POINTS 

8240 

WEIGHTING  FACTORS 

8240 

REFERENCE  SEQUENCES 

10000 

RAW  DATA  FOR  FFT 

1024 

TOTAL 

27504  WORDS 

MEMORY  REMAINING  FOR 
PROGRAM  EXECUTION 

5264  WORDS 

Figure  10.  Memory  liequirements 


an  alternative  will  be  discussed  later.  Once  again,  the  estimates  provided  in  this 
figure  are  not  exact.  Rather,  they  are  estimates  which  show  that  memory  limita- 
tions appear  to  be  more  severe  than  processing  time  limitations  in  implementing 
this  form  of  TUSAR  processing.  Stress  should  then  be  placed  on  programs  which 
economize  memory  should  actual  requirements  dictate. 

Figure  1 1 takes  a look  at  a block  diagram  representation  of  memory 
allocations. 

It  is  important  to  remember  that  evaluated  processing  time  and  storage 
requirements  apply  to  200  moving  target  doppler/ doppler-rate  combinations  in  a 
single  range  "bin.  " At  a distance  of  400  m,  we  will  be  confined  to  observing 
objects  within  an  18  M range  resolution  cell.  For  other  range  resolution  cells, 
one  must  add  another  dimension  to  the  time  and  memory  requirement  calculations 
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Figure  11.  Data  Flow  and  Memory 
Hequirements 


27  Hilt  I Vt  H 


or  reduce  some  of  the  200  dopplers/doppler-rate  combinations/  in  favor  of  addi- 
tional range  coverage.  For  the  truck  case,  it  will  generally  he  necessary  to 
consider  only  one  range  interval  at  a time. 


12.  SHOKTK.lt  MKMOIO 

One  unique  method  for  cutting  down  on  memory  requirements  is  to  use  a 
triangular  rather  than  rectangular  array  of  data  storage.  Once  the  size  of  the 
F FT  has  been  established,  there  are  a fixed  number  of  frequency  bins  and  time 
periods  to  be  considered.  In  the  previous  memory  calculations,  a rectangular 
array  of  fixed  size  was  used  to  compute  the  memory  size,  and  was  valued  at 
80  x 103  words.  This  rectangular  array  of  data  assumes  a need  to  keep  every 
output  from  all  integration  time  intervals  for  all  calculations,  which  is  not  neces- 
sarily the  case.  Because  of  the  linear  decrease  in  target  doppler  frequency  with 
time,  the  later  in  time  an  FFT  output  is  required  for  a particular  coherent  sum- 
mation, the  less  likely  will  be  the  requirement  for  a high-frequency  output.  It 
becomes  clear  that  reduced  memory  results  from  saving  only  those  low-frequency 
FFT  outputs  likely  to  be  used.  If,  after  FFT  processing,  new  FFT  output  data  is 
stored  onto  the  high-frequency  data  points  of  previous  FFTs,  a triangular  matrix 
results  that  requires  25%  to  50";.  less  memory.  The  reduced  memory  formation 
requires  additional  memory  and  time  for  keeping  track  of  all  the  data  locations 
presently  in  the  memory  structure. 

A highly  simplified  version  of  the  triangular  storage  format  is  displayed  in 
Figure  12.  Within  ench  triangular  memory  array,  a column  represents  a specific 
FFT  frequency  output.  Numbers  represent  the  time  sequence  during  which  the 


17 


— m 


10* 

4 

mad 

(0* 

med 

h.gh 

high 

3 

4 

2 

3 

4 

0 

® 

® 

m 

4 

0 

L 

0 

0 

6 

6 

6 

5 

5 

0 

u 

0 

7 

7 

L 

6 

0 

® 

0 

0 

®i 

4 

3 

4 

® 

® 

10 

5 

3 

3 

0 

0 

b. 

7 

7 

6 

6 

0 

5 

j 

® 

7 

1 

8 

1 

7 

® 

® 

9 

9 

8 

9 

® 

t. 


Figure  12.  Triangular  Data 
Structure  (Simplified).  Circled 
numbers  represent  outputs 
summed  that  time  period, 
columns  are  frequencies, 
numbers  are  FFT  number 
(highest  most  recent). 


FFT  outputs  are  generated  and  stored.  Circled  numbers  represent  data  which  is 
combined  to  produce  the  detection  statistic  and  hence  memory  locations  available 
for  new  data.  In  progressing  from  (a)  to  (f)  in  Figure  12,  it  is  clear  how  new  FFT 
data  replaces  data  which  has  just  been  used.  Only  the  low-frequency  outputs, 
which  will  be  used  as  the  FFT  grows  old  and  the  time  axis  shifts,  are  saved  the 
whole  time.  The  higher  the  frequency,  the  sooner  it  is  used  and  discarded.  Of 
course,  this  procedure  will  work  best  if  the  FFT  size  is  predetermined,  allowing 
the  programmer  to  be  sure  that  certain  outputs  can,  indeed,  be  overwritten. 

In  greater  detail,  here  is  how  the  data  structure  of  Figure  12  works.  The 
data  from  FFT  is  initially  set  along  diagonals,  and  counters  for  each  frequency 
column  are  set  to  the  base  of  the  triangle.  Each  time  period,  circled  data  points 
are  used.  These  memory  locations  become  available  to  store  the  next  FFT  outputs. 

Note  in  Figure  12  that  the  column  counters  only  need  to  count  up  to  the  top  of 
each  column,  reset,  and  begin  counting  up  again.  Thus,  one  needs  only  as  many 
counters  as  one  has  frequencies. 

The  properties  of  the  triangular  memory  structure  as  assumed  target  param- 
eters change  are  of  interest.  The  closer  the  target  slope  m is  to  -1,  the  more 
triangular  is  the  structure.  As  the  slope  shallows,  the  data  structure  appears 
more  like  a cornered  rectangle  (Figure  13).  This  is  due  to  the  fact  that  at  shallow 
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Figure  13.  Data  Structures  for 
Different  Slopes 
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slopes,  a given  frequency  will  be  used  a few  times,  rather  than  once.  At  steeper 
slopes,  several  frequencies  may  be  thrown  away  at  once,  so  an  even  smaller  tri- 
angle is  required  in  terms  of  relative  dimensions,  but  the  number  of  data  points 
included  increases,  so  that  a tradeoff  does,  indeed,  occur. 


13.  CONCLUSION 

It  is  clear  that  the  design  of  an  efficient  TUSAR  processor  is  not  a simple  task. 
Care  must  be  taken  to  trade  off  integration  time,  memory  and  computation  load  in 
selecting  the  appropriate  FFT  size.  Interface  between  the  radar  video  data  and 
processor,  and  between  the  summation  output  and  display,  must  be  considered. 

One  must  also  consider  the  magnitude  of  the  problem:  a postulated  target  can  have 
any  velocity  vector  magnitude  and  direction,  and  the  radar  platform  can  look  at 
any  of  a number  of  ranges  for  a variety  of  platform  velocities,  yielding  a four- 
dimensional array  of  reference  data. 

Despite  the  obvious  complexities,  the  project  is  by  no  means  hopeless.  Time 
and  memory  restrictions  do  not  severely  limit  laboratory’  project  scope  or 
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capabilities.  Programming  parameters  are  not  highly  sensitive  to  target  and 
radar  parameter  changes,  particularly  if  a triangular  memory  allocation  system 
is  used. 

In  the  final  analysis,  an  operational  model  will  be  an  order  of  magnitude  more 
complex  than  the  laboratory  model  of  the  radar.  As  few  pilots  are  willing  to  fly 
at  a constant  speed  over  a straight  course  for  any  length  of  time  in  a combat 
environment,  larger  memory  is  required  along  with  the  use  of  plug-in  PROM 
cards  to  accommodate  different  aircraft  speeds.  The  use  of  the  TUSAR  technique 
on  board  a C-130  or  similar  aircraft  is  realistic,  and  apparently  lies  within  our 
progressing  digital  processing  capabilities. 
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Appendix  A 


Listing  and  Sample  Outputs 


This  Appendix  contains  a listing,  and  sample  outputs  for  the  computer  pro- 
gram described  in  Section  8. 
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«input  .OUTPUT  .T.Ptfc)  ;SJfM 
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r.wn,$ 
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TINC»?*PX*1/0FIM 

RESOLVE  THE  ANGLES  TO  RAOIANS 

- RTA»PT*RTA/180 
VTI«Pl*VTA/teO 

FINO  THE  RECTANGULAR  CCORPINATES  OF  PELT  A V BAR  ANO 
THETA  SUB  0 OEL  V. 

AOELV*VR-VTM*COS (VTA) 

B0€LV«-VTH*SIM( VTA) 

- ANG0V*RTA-ATAN2(80flV,ADElV>  — 

fimo  the  target  parameters 

MO«V*PI*FR*SQPT( AOELV**240CEL V**2> *COS (ANGOV) / ( 3M0**BI 
WD0TT*-4*PI*EP* ( A0ELV**?*80ELV*#?)*SIN (ANGOV) ##2/(RTH#3#lO##8) 

ENO  OF  PARAMETER  EINOING.  ROM  BEGIN  SEARCH  ROUTINE.* 


()n*M20 
flOO 130 
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0*0  310 
000360 
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000400 

000410 

00"42Q 
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000440 

0004*0 

000460 

000470 

0004*0 

mm 
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000520 
000530 
000540 
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000570 

mm 
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mm 
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N*  1 f IZ»0 

MFrWO-0,  ?*MCM  £ LF«INT  ( ME/ DEL  Ml 

IF  IMF.it. 0)  l*«LF-l  

000670 

CALCULATE  ENTERING  ANC  EXITING  FREQUENCIES  FOP  GIVEN  TINE  INTERVAL0 0 06 6 0 
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iBiKCKVlhUift^tk  °ono0oH°o 

MSTOP*MOOTT*TlNC*N*MO  000730 

000740 

FINO  INTEGER  VALUE  OF  EXIT  FREQUENCY  INTERVAL  J0075jj 
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000960 
000860 
000*6 0 
000870 
000890 
000690 

0 09  0 0 
000)10 
000920 
000930 
go  0Q40 
000950 
000960 
000970 
porq*n 
000990 
001 000 
001010 

mm 

mm 

001060 

001  07  0 
001080 
0 0 1 0 ° 0 
001100 
001110 
001120 
001 1 *0 
001160 
001150 
001160 
001170 
cot  190 
001190 
001200 
001210 

15155* 

*)001 760 
001250 

0 0 1 ? F 0 
nnt  ?70 
001280 

*515?* 

MUM 

001330 

001  160 

ssiks 

881558 

mm 

001610 

001620 

001630 

001660 

001650 

001660 

001670 

001690 

0016°0 

001600 

001610 

001520 

001630 

nm* 

001560 
001570 
001680 
PM59D 
001600 
001610 
001620 
’ 0 


23 


TARGET  is  AT  RANGE  50000. C METERS*  bearing  90.0  OEGREES  PROM  FLIGHT  PATH. 

TARGET  VELOCITY  IS  3.66  METERS/SEC  AT  RELATIVE  BEARING  45.0  OEGREES  (RIGHT  MING  POSITIVE). 
LOP  RATE  IS  50.0  PERCENT.  

RAOAR  IS  SET  UP  AS  FOLLOMSN 

VELOCITV\  154.33  METERS/SEC. 

PRF\  250  HERTZ. 

CARRIERS  • 100E»09  HERTZ. 

EACH  FFT  IS  512  POINTS  LONG.  

WO*  -.02  TIMES  HCM. 

MOOTT  * .96  TIMES  WOOTCm. 


ME 

REQUIRE  39  OUTPUTS 

AS  FCLLOMSN 

POSITION  FFT 

OUTPUT  • 

0 

508  . 

508 

1 

509 

5 09 

2 

1019 

507 

3 

1529 

505 

4 

1530 

506 

5 

2040 

504 

6 

- ---  2551  _ 

503 

7 

3062 

502 

0 

3572 

500 

9 

3573 

501 

10 

4083 

499 

11 

4594 

498 

12 

5104 

496 

13 

5105 

497 

14 

5615 

4 95 

15 

6126 

494 

16 

6637 

493 

17 

7147 

491 

10 

7148 

49  2 

19 

7658 

490 

20 

8169 

489 

21 

8679 

487 

22 

8680 

488 

23 

9190 

486 

24 

9701 

485 

25 

10212 

4 84 

26 

10722 

482 

27 

10  723 

483 

20 

11233 

481 

29 

11744 

480 

30 

12254 

478 

31 

12255 

479 

32 

12765 

477 

33 

13276 

476 

34 

13787 

475 

35 

14297 

473 

36 

14298 

474 

37 

14808 

472 

30 

15319 

471 

31  CHOPS  MAOE. 

39  OUTPUTS  LEFT. 

30  TIME  PERIODS  ELAPSEO. 
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t I 


n 


„B  js  best  «b*wh 

moll  OOH  EUB»1SBS1  TO  U>C  " 


TARGET  IS  AT  RANGE  5 0000.0  HETERS,  B EARING  99.0  OEGRCES  FROM  FLIGHT  PATH. 

TARGET  VELOCITY  IS  3 . 86  HETERS/SEC  AT  RELATIVE  REARING  0.0  OEGRECS  (RIGHT  MING  POSITIVE), 
LOP  RATE  IS  90.0  PERCENT. 


RAOAR  IS  SET  UP  AS  FOLLOMSV 

VELOCITTN  159.33  HETERS/SEC. 
PRFV  ?50  HERTZ. 

CARRIERS  • 1 00E ♦ 09  HFRT 7 • 

EACH  FFT  IS  51?  POINTS  LONG* . 
MO*  -.00  TINES  MCH. 

MOOTT  • .95  TIHES  MOOTCH. 


ME  REQUIRE  39  OUTPUTS 
POSITION  FFT 

51?  

10?? 

10Z3 
1533 
?099 
?555 

3065 

3066 
3576 
**067 
(.596 
5106 

5109  

5619 
6130 
6691 
7151 
715? 

766? 

0173 
0661 
0609 
9199 
9 709 

10  ?16 

10  7?6 
10  7?7 
11237 
11796 
1??59 

12769  

12770 
13260 
13791 
1930? 

1901? 

19013 
15  3?3 
15639 


AS  FCLL09SN 
OUTPUT  0 
51? 

510 

511 
509 

906 
507 

505 

506 
509' 

503 

50? 

500 

501 
999 

996 

997 

995 

996 
999 
99  3 
991 
99? 

9 90 

969 
966 
9 66 

907 
965 
909 
963 
901 
96? 

960 

979 

970 
9 76 
977 
975 
979 


3?  CHOPS  HAOE. 

39  OUTPUTS  LEFT. 

31  TINE  PERIOOS  ELAPSED. 
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target  is  at  range  socoo.o  meters,  rearing  so.#  degrees  from  flight  path. 

TARGET  VELOCITY  IS  3. Oft  METERS/SEC  AT  RELATIVE  BEARING  1J5.0  DEGREES  (RIGHT  MING  POSITIVE). 
LOP  RATE  IS  $0.0  PERCENT. 

RAO AP  IS  SET  UP  A$  FOLLCMSV 

VELOCITY\  ISA. J1  NETERS/SFC. 

PRE  \ 3S0  HERT7, 

CAPRIERV  . 1 00E  .09  MERT7. 

EACH  FFT  IS  $17  POINTS  LONG. 

NO*  -.07  TINES  M CM. 

NOOTT  ■ 1.0$  TINES  MOOT CM. 

HE  REQUIRE  39  OUTPUTS  AS  FOLLOMSN 


POSITION 

FFT  OUTPUT  • 

n 

508 

506 

1 

509 

50  9 

2 

1019 

607 

1529 

605 

h 

1530 

506 

9 

2060 

506 

2550 

502 

7 

2 551 

503 

6 

3061 

501 

9 

357? 

500 

10 

608? 

696 

i 11 

6083 

699 

1 " 

6693 

69  7 

13 

5103 

695 

1% 

5106 

696 

**  19 

5616 

696 

16 

6125 

693 

i if 

6635 

691 

18 

6636 

69? 

19 

7166 

690 

20 

7657 

689 

i 2 1 

8167 

687 

» 27 

8166 

686 

23 

OhlO 

666 

2k 

9186 

666 

} 29 

9189 

685 

26 

9699 

663 

27 

10210 

66  2 

20 

10  720 

680 

29 

10  721 

661 

30 

11231 

679 

31 

11762 

676 

32 

12252 

676 

33 

12253 

677 

36 

12763 

675 

35 

13273 

67  3 

36 

13276 

676 

37 

13786 

672 

38 

16295 

671 

29  CHOPS  H60E. 

39  OUTPUTS  LFFT. 

70  TINE  PERIODS  ELAPSEO. 


I 


k 


moil  CQi’¥  FURNISHED  xo  - 


t% 

A U* 

1 76 

lAROrT 

tS  AT  RANGE 

'0000.0 

targct 

VtlOCHY  IS 

3.06  MfTl 

io.  *»rc  is  i».«  cr.ctNi. 

RAOAR  IS  SF I UR  AS  l 

FOIL  QMS \ 

VC LOCXTtX  ISt, 3 T 

hctcrs/: 

MF\ 

?tO  NCRT7 

.tOOF.M  HFKII. 

f ACM 

FFT  I«  70tS 

POINTS  1 

MO- 

i •• 

• 7 IXMfS  MCM 

MOOT!  • 

.st  Unrs 

nootcn. 

Mf 

trOUXRF  16t  OUTPUTS  AS  F| 

position 

FFT  OUTI 

• 

?0lt 

7tlt 

1 

701% 

70  It 

7 

701 1 

7016 

1 

TUT 

T017 

A 

7 0 1 A 

7010 

% 

701% 

?on 

• 

7370 

7070 

7 

7071 

7071 

6 

7077 

7077 

% 

7073 

7073 

10 

7 0?t 

?0?t 

II 

70? % 

7070 

1? 

7076 

7076 

I) 

7077 

707? 

1% 

7070 

7070 

IS 

7 0.*% 

707% 

It 

7030 

70  TO 

IT 

7031 

7031 

It 

7037 

703? 

IS 

70  IS 

7033 

tl 

7 0 3t 

70  3t 

ft 

tttl 

19%  1 

If 

tot? 

l%%t 

tl 

tots 

l%%t 

l % 

tott 

1996 

7% 

tot% 

l%%? 

1 1 

tott 

1%%9 

tr 

tot? 

1S99 

ti 

toto 

7000 

2% 

tott 

7001 

St 

toto 

7007 

Jl 

toil 

7001 

It 

t0%7 

?00t 

33 

tots 

7009 

St 

tott 

7006 

St 

tott 

700? 

St 

tott 

7000 

IT 

tot?  — 

700% 

St 

tot* 

7010 

si 

tott 

7011 

to 

toto 

701? 

tl 

totl 

7013 

1 1 

6066 

1 %7? 

ts 

tott 

I%7% 

tt 

(.070 

l%?t 

t* 

6071 

1979 

tt 

60  77 

1976 

t? 

607S 

1%77 

tt 

60  7 1 

1%76 

tt 

t«Tt 

%V% 

► ATM. 

CRX&MT  MING  POSXTIVf) , 
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TARGET  IS  AT  RANGE  50000.  0 NET  ERS»  BEARING  50.0  OEGREES  FROM  FLIGHT  PATH* 

TARGET  VEIOCITV  IS  3.86  NETfRS/SEC  AT  RELATIVE  SCARING  UO.O  OEGREES  (RIGHT  MING  POSITIVE), 
LOP  RATE  IS  50,0  PERCENT.  

RAOAR  IS  SET  UP  AS  FOltOMSX 

VELOCITTN  154.33  HETf RS/Sf C. 

PRF\  750  HFRTZ • 

CARRIE R\  • 100E>09  HERTZ. 

EACH  FFT  IS  *56  POINTS  LONG. 

MO»  *00  TINES  MCH. 

MOOT!  ■ 1.05  TINES  WOOTCN. 

ME  REQUIRE  54  OUTPUTS  AS  FOltOMSX 
POSITION  FFT  OUTPUT  • 

0 *56  756  . 


l 

517 

756 

7 

1073 

755 

3 

! 779 

?«S 

6 

1790 

754 

5 

201.6 

754 

6 

7557 

753 

t 

?St3 

293 

a 

3068 

75? 

9 

3 069 

753 

10 

3374 

75? 

11 

3500 

75? 

1? 

3835 

761 

13 

3836 

75? 

14 

4 091 

751 

15 

4347 

751 

16 

460? 

750 

17 

4603 

751 

IS 

4658 

750 

19 

5114 

750 

?0 

5369 

749 

71 

5370 

?60 

22 

5625 

749 

23 

6136 

740 

2k 

6 39? 

748 

2 5 

6903 

747 

To 

7 169 

747 

27 

7670 

290 

26 

7976 

746 

2 9 

8181 

745 

30 

816? 

746 

31 

0937 

745 

32 

669T 

745 

33 

6948 

744 

34 

8949 

745 

35 

9704 

244 

36 

9460 

744 

37 

9715 

743 

30 

9 716 

744 

39 

9971 

74  3 

40 

10227 

743 

41 

10982 

24? 

47 

10403 

743 

43 

10  730 

74? 

44 

11249 

741 

45 

11505 

741 

46 

12010 

240 

47 

12272 

740 

40 

12783 

739 

739 

26 


M 
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Printed  by 

United  States  Air  Perce 
Hentcem  APB,  Mess.  01731 


i 


